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Abstract. Three-nucleon interactions are a frontier in understanding and predicting the structure of 
strongly-interacting matter in laboratory nuclei and in the cosmos. We present results and discuss 
the status of first calculations with microscopic three-nucleon interactions beyond light nuclei. This 
coherent effort is possible due to advances based on effective field theory and renormalization group 
methods in nuclear physics. 



The physics of strong interactions extends over extremes in density, neutron-to-proton 
imbalance and temperature (see Fig. 1). Understanding and predicting the properties 
of these fascinating forms of matter requires progress on fundamental problems in the 
theory of nuclear forces and in many-body physics. In this talk, we highlight the key role 
of three-nucleon (3N) interactions for nuclear structure. 

Nuclear interactions depend on a resolution scale, which we denote by a generic 
momentum cutoff A, and the Hamiltonian is always given by an effective theory for 
nucleon-nucleon (NN) and corresponding many-nucleon interactions (with correspond- 
ing effective operators) [1, 2, 3]: 

H(A) =T + Vnn(A) + V 3 n(A) +V 4 n(A) + . . . . (1) 

This scale dependence is similar to the scale dependence of parton distribution functions, 
and shows that the effect of 3N interactions depends on this scale and the theoretical 
convention that enters all parts of the Hamiltonian. 

At very low momenta Q <C m n ~ 140 MeV, the details of pion exchanges are not re- 
solved and nuclear forces can be systematically expanded in contact interactions and 
their derivatives [1]. The corresponding pionless effective field theory (EFT) is ex- 
tremely successful for capturing universal large scattering-length physics (with improve- 
ments by including effective range and higher-order terms) in loosely-bound or halo nu- 
clei, reactions at astrophysical energies and in dilute neutron matter (see e.g., Ref. [4]). 

For most nuclei, the typical Fermi momenta are Q ~ m n and therefore pion exchanges 
have to be included explicitly. In chiral EFT, nuclear interactions are then organized 
in an expansion in powers of Q/A^, where At, denotes the breakdown scale, roughly 
Ab ~ mp [1,2]. The great advantage is that up to next-to-next-to-next-to-leading order 
(N 3 LO) only two new 3N couplings enter, since parts of the 3N force can be consistently 
determined by tzH and NN couplings. The two 3N couplings, as well as the short-ranged 
NN couplings, depend on the resolution scale A and for each A can be fit to data. 

Using the renormalization group (RG), we can evolve an initial potential to lower 
resolution by integrating out high momenta through discretized RG flow equations 
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FIGURE 1. From the QCD vacuum to laboratory nuclei to neutron stars and supernovae. 
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FIGURE 2. Evolution of the 3 Si partial wave with a smooth n eX p = 8 regulator at cutoffs A = 



5.0,4.0,3.0,2.0 and 1.5fm 1 (for details see Ref. [6]). The initial potentials are Argonne Vj8 [11] (top) 
and a N 3 LO chiral EFT potential from Ref. [12] (bottom). The color scale ranges from —0.5 to 0.5fm. 



or equivalent Lee-Suzuki transformations [3, 5]. In the last two years, these methods 
have been refined to employ smooth regulators [6] and similarity renormalization group 
(SRG) transformations [7, 8], which both have technical advantages in oscillator spaces. 

Changing the cutoff leaves observables unchanged by construction, but shifts contri- 
butions between the potential and the sums over intermediate states in loop integrals. 
Since the sums are restricted by the intrinsic resolution A, these shifts can weaken or 
largely eliminate sources of nonperturbative behavior such as strong short-range repul- 
sion and short-range tensor forces [9, 10]. As shown in Fig. 2, the evolution leads to 
low-momentum interactions, generally known as "Viowfc", which become universal for 
A < 2fm~ 1 and have weak off-diagonal coupling. Since the RG preserves the long-range 
parts of nuclear interactions [6], this leads to a weak renormalization of long-range op- 
erators, for example for the deuteron rms radius in Fig. 3. 

The evolution of chiral EFT interactions to lower cutoffs is beneficial for two reasons. 
First, the RG generates all higher-order short-range contact interactions so that observ- 
ables are exactly reproduced and the theoretical uncertainty remains at the level of the 
truncation error in the initial potential. We illustrate this in Fig. 4 by comparing the flow 
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FIGURE 3. The deuteron rms radius rj calculated with the bare operator as a function of the cutoff 
(left). The radius is only weakly dependent on the cutoff for A > 1 .5 fm~ 1 , although the intrinsic resolution 
in the S-state deuteron wave function changes substantially with A (right). For details see Ref. [6]. 
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FIGURE 4. Flow of Vi ow i c (k / — Q,k = 0;A) compared to the corresponding momentum-independent 
contact interaction Co (A) at LO and NLO, where this coupling is determined entirely from RG invariance 
and fits to the scattering length a s (at LO) plus effective range r e (at NLO). 

of Viowkft = 0,fc = 0;A) with the corresponding momentum-independent contact in- 
teraction Co (A) in subsequent orders of pionless EFT. Second, lower resolutions can be 
represented efficiently in oscillator bases, and therefore lead to direct convergence in nu- 
clear structure calculations [6]. This is demonstrated by the very promising convergence 
for N max ~ 10 in NCSM calculations with SRG interactions [13] shown in Fig. 5. 

Viowjfc(A) defines a class of NN interactions with cutoff-independent low-energy NN 
observables. Consequently, any cutoff variation of observables estimates the truncation 
errors due to neglected many -body interactions in H(A) of Eq. (1). Figure 6 shows 
that this cutoff dependence explains the empirical Tjon line, that 3N interactions are 
required by renormalization and needed to break off the line to accurately describe 
the experimental 3 H and 4 He binding energies [14]. NN-only results also lead to Tjon 
lines in medium-mass nuclei, where results truncated in oscillator shells as well as for 
different HQ. lie approximately on the same lines [15]. 

Three-nucleon interactions are crucial for many-nucleon systems. In addition to the 
impact on binding energies, they play a central role for spin-orbit and spin dependences, 
for isospin dependences of neutron- and proton-rich systems, and they drive the density 
dependence of nucleonic matter. The latter are pivotal for extrapolations to the extremes 
of astrophysics. Three-nucleon interactions are also a frontier in few-body scattering (for 
a critical discussion, see e.g., Ref. [16]). Since 3N contributions are amplified in nuclei, 
it may be necessary for controlled predictions to constrain 3N couplings with few- and 
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FIGURE 5. Ground-state energy of 6 Li as a function of the oscillator parameter hD. at four different 
values of SRG resolution A = 3,2, 1.5 and 1 fm _1 . The NCSM results clearly show improved convergence 
with the maximum number of oscillator quanta N max for lower cutoffs. For details see Ref. [13]. 
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FIGURE 6. Correlation of the 3 H and 4 He binding energies. The cutoff dependence of the exact NN- 
only results with Vi OW A:(A) explains the empirical (solid) Tjon line. For details see Ref. [14]. 
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FIGURE 7. Leading (N 2 LO) 3N interaction in chiral EFT without explicit Deltas. 



many -body data. Therefore a coherent 3N effort is needed with theoretical uncertainties. 

In chiral EFT without explicit Deltas, 3N interactions start at N 2 LO or order 
(<2/Ab) 3 [1, 2] and their contributions are given diagrammatically in Fig. 7. The 
low-energy constants for the long-range parts relate 7rN, NN and 3N interactions, and 
the determination from 7rN scattering is, within errors, consistent with the extraction 
from peripheral NN waves. The present constraints are c\ = — 0.9l[j'g, C3 = — 

and C4 = 3.5^Q 2 (all in GeV -1 ) [17]. In particular, C3 and C4 are important for nuclear 
structure and have large uncertainties (see Fig. 1 1 for the impact on nucleonic matter). 

For lower cutoffs, we take the corresponding 3N interactions V3n(A) from chiral 
EFT by fitting the D- and E-term couplings to the 3 H and 4 He binding energies for a 
range of cutoffs [14]. Since chiral EFT is a complete basis, this gives the 3N force up to 
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FIGURE 8. The J H and 4 He radii are approx. cutoff independent with NN and 3N interactions [18]. 
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FIGURE 9. Relative contributions | AE /£ | to the binding energy of 4 He at the CCSD level from V lowjfe as 
well as normal-ordered 0-, 1-, 2-body and residual 3-body parts of the 3N interaction (left). Convergence 
of the corresponding CCSD(T) results with the number of oscillator shells (right). The extrapolated 
binding energy — 28.23 MeV agrees well with the exact result — 28.20(5)MeV. For details see Ref. [19]. 



truncation errors. We have found that the resulting 3N interactions become perturbative 
for A < 2fm~ 1 , while they are nonperturbative for larger cutoffs, and that the size of 3N 
expectation values is natural ~ (2/Ab) 3 (V / i OW yt) [14]. 

As for the Tjon line, cutoff variation can provide lower bounds for theoretical uncer- 
tainties due to neglected many-body interactions or an incomplete many-body treatment. 
As shown in Fig. 8, this is a powerful and practical tool, which is important for extrap- 
olations, in particular for matrix elements needed in fundamental symmetry tests (e.g., 
double-beta decay and isospin-violating corrections for super-allowed beta decay). 

Coupled-cluster (CC) theory combined with rapid convergence for low-momentum 
interactions pushes the limits of accurate calculations to medium-mass nuclei and sets 
new benchmarks for 16 and 40 Ca [19]. In Fig. 9, first CC results with 3N forces show 
that low-momentum 3N interactions are accurately treated as effective 0-, 1- and 2-body 
terms, and that residual 3N interactions can be neglected [19]. This is very promising 
and supports the idea that phenomenological monopole shifts in shell model interactions 
are due to 3N contributions [20], see also the talk by T. Otsuka and Ref. [21]. This links 
understanding the shell model and the drip lines to 3N forces. First investigations in this 
direction are presented in Fig. 10 [22]. We speculate that the T = monopoles with 
large cutoff dependences receive cutoff-dependent contributions mainly from second- 
order NN-3N terms (which are expected to be attractive) and that the T = 1 monopoles 
receive cutoff-independent contributions from the 3N c ; terms (which are repulsive in 
nuclear matter). Systematic investigations of 3N forces and valence-shell interactions 
are in progress [22]. 

Finally, there is the possibility of perturbative nuclear matter with low-momentum 




FIGURE 10. Monopole interactions in the pf shell obtained from Vi 0W £ plus second-order many -body 
contributions for a range of cutoffs, compared to phenomenological matrix elements [22]. 
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FIGURE 11. Nuclear (left) and neutron matter (right) from low-momentum NN and 3N interactions 
based on Hartree-Fock plus dominant second-order contributions for various cutoffs. For details see 
Ref. [9, 23]. 



NN and 3N interactions [9]. As shown in Fig. 11, with second-order contributions the 
cutoff dependence is very weak at low densities. We have found that 3N forces drive 
saturation [9], and at present the uncertainties in the c ; overwhelm other uncertainties 
(see the results for neutron matter) [23]. Lower values for C3 (as also expected from 
N 3 LO 3N terms) add w — 2MeV at saturation density in neutron matter, and this also 
corresponds to the value from Delta resonance saturation (which explains the agreement 
with the results of Akmal et al.). These findings will provide guidance for constructing 
a universal density functional for nuclei. Towards denser matter, 4N force contributions 
of E/A ~ 1 MeV would not be unreasonable. 

This is an exciting era: with advances on many fronts, a coherent effort to understand 
and predict nuclear systems based on effective field theory and renormalization group 
interactions, where 3N forces are a frontier. 
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